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ABSTRACT 

We report the results of a spectroscopic search for Lyman-a (Lya) emission from gamma-ray 
burst (GRB) host galaxies. Based on a well-defined parent sample (the TOUGH sample) of 69 
X-ray selected Swift GRBs, we have targeted the hosts of a subsample of 20 GRBs known from 
afterglow spectroscopy to be in the redshift range z — 1.8-4.5. We have obtained spectroscopy 
using the FORS1 instrument at the ESO Very Large Telescope to search for the presence of Lya 
emission from the host galaxies. We detect Lya emission from 7 out of the 20 hosts, with the 
typical limiting 3c line flux being 8 x 10 -18 erg cm -2 s" 1 , corresponding to a Lya luminosity of 
6 x 10 41 ergs" 1 at z = 3. The Lya luminosities for the 7 hosts in which we detect Lya emission 
are in the range (0.6-2.3) x 10 42 ergs -1 , corresponding to star- formation rates of 0.6-2.1 M Q yr 
(not corrected for extinction). The rest-frame Lya equivalent widths (EWs) for the 7 hosts are 
in the range 9-40 A. For 6 of the 13 hosts for which Lya is not detected we place fairly strong 
3cr upper limits on the EW (< 20 A), while for others the EW is either unconstrained or has 
a less constraining upper limit. We find that the distribution of Lya EWs is inconsistent with 
being drawn from the Lya EW distribution of bright Lyman break galaxies (LBGs) at the 98.3% 
level, in the sense that the TOUGH hosts on average have larger EWs than bright LBGs. We 
can exclude an early indication, based on a smaller, heterogeneous sample of -pie-Swift GRB 
hosts, that all GRB hosts are Lya emitters. We find that the TOUGH hosts on average have 
lower EWs than the pie-Swift GRB hosts, but the two samples are only inconsistent at the 92% 
level. The velocity centroid of the Lya line (where detected) is redshifted by 200-700 km s" 1 with 
respect to the systemic velocity (taken to be the afterglow redshift), similar to what is seen for 
LBGs, possibly indicating star-formation driven outflows from the host galaxies. There seems to 
be a trend between the Lya EW and the optical to X-ray spectral index of the afterglow (Pox), 
hinting that dust plays a role in the observed strength and even presence of Lya emission. 



Subject headings: dust, extinction — galaxies: fundamental parameters 
galaxies: star formation — gamma-ray burst: general — surveys 

1. Introduction 



galaxies: high-redshift 



Based on observations collected at the European 
Southern Observatory, Chile, as part of program 177. A- 
0591. 
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Due to their potential brightness at wavelengths 
ranging from radio to gamma-rays, gamma-ray 
bursts (GRBs) and their afterglows can be used 
as powerful astrophysical probes. They are mo- 
mentarily bright enough to be observed any- 
where in the Universe, even if located in the 
most dusty environments or at the highest red- 
shifts (e.g. IWiiers et al. 1998), but later fade 
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away and allow a detailed study of their envi- 
ronment. It is now established that long-duration 
GRBs are associated with core collapse super- 
novae (e.g.JStanek et aLll2003t iHiorth et al.ll2005 



Wooslev fe Blooml 120061: IHiorth fe Blooml 1201 

and hence star-formation. GRBs thus can be used 
to probe the star-formation density over most of 
the cosmic history from the formation of the first 
stars to the present. The most spect acular exam- 
ple o f this is GRB 090423 at z = 8.2 dTanvir et al 



20091 : Salvaterra et al. 20091) representing a look- 
back time of more than 95% of the time since the 
Big Bang. 

Lya-emitting galaxies are in a state of active 
star-formation and most likely contain little or 
no dust, since Lya photons have a much higher 
probability than other UV photons of be ing ab- 
sorbed by dust due to resonant scattering (lAdamsl 



19931) 



Il972t ICharlot fe Fall 11993k IValls-Gabauc 
Lya emission is hence sensitive to both the star- 
formation rate and the dust content of GRB host 
galaxies. In addition, geometrical effects and the 
kinematical state of the interstellar medium seem 
to be important for the esc a pe of Lya photons 



g., Giavalisco et al. 19961: Haves et al. 20051: 



Verhamme et all 120061 iLaurseil et al.ll2009h . 



The first studies of Lya emission from (pre- 
Swift) GRB host galaxies indicated that Lya emis- 
sion seemed to be ubiquitous, with 5 detections 
out of 5 possible (|Fvnbo et al.ll2003[ f T l. This would 
be intriguing as only about 25% of Lyman-break 
selected galaxies (LBGs) at similar redshifts have 
Lya emission with res t-frame equivalent width 
(EW) larger than 20 A (|Shaplev et al.ll2003l): this 



is als o the case at higher redshifts (IDouglas et al 



120101 but see also IStark et~al"1l20lC Another 6 
GRB host Lya emitters (all from Swift GRBs) 
have been reported in the literature since therQ 
(excluding the hosts reported in this worljf]), but 



1 GRB971214 llKulkarni et al.l Il998h. GRB 00092 6 
l|Fvnbo et ail l2002h. GRB 0112 11 l|Fvnbo et al.l l2003h . 
GRB 021004 llMeiler et al.l |2002|. and r eferences therein), 
and GRB 030323 [|Vreeswiik et al.ll2004f ). 

2 GRB 060714 l lJakobsson et ail l2006bl ). GRB 060926 
jFvnbo et ail 120091) GRB061222A l|Perlev et ail 
120091). GRB 07011 jFvnbo et ail l2009t). GRB 071031 
llFvnbo et ail l2009f ). and GRB 090205 jp'Avanzo et ail 
2010). 

3 This work reports the detection of Lya emission from 7 host 
galaxies, of which one (GRB 070110) was already identified 
as a Lya emitter in the literature. 



there still has not been a systematic examination 
of the frequency of Lya emitters among GRB host 
galaxies. This is the aim of the present work. 

Given the effect of dust on Lya photons, pos- 
sible explanations for an excess of Lya emitters 



amon g GRB host galaxies include ([Fvnbo et al 



2003): (i) a preference for GRB progenitors to 
be metal poor as expected in the collapsar m odel 
dWooslev fc Hegerll2006l:[Yoon fc Langerll2005l see 



also lNiino et al.ll2009l ): (ii) an optical afterglow se- 
lection bias against dusty hosts; (iii) a higher frac- 
tion of Lya emitters at the faint end of the high- 
redshift luminosity function, where most GRB 
hosts arc found; (iv) small-number statistics. Us- 
ing a well selected and more complete Swift sample 
we shall here address these issues. 

The paper is structured in the following way. 
The parent sample (TOUGH), the target selec- 
tion, the spectroscopic observations and the data 
reduction are described in £]3J The Lya detections 
and upper limits are presented in £13.11 The veloc- 
ity offset of the Lya emission with respect to the 
systemic velocity as given by the afterglow red- 
shift is discussed in £13.21 A comparison of the 
Lya fluxes from afterglow and host spectra is done 
in £13.31 In £j4] we discuss the results, including 
how they relate to LBGs and to pre- Swift stud- 
ies, and how the observed Lya emission is related 
to the afterglow broad-band spectral index /?ox, 
and we summarize our findings. Finally, the Ap- 
pendix presents observations targeting the hosts of 
3 GRBs that are not part of the complete, well de- 
fined TOUGH sample discussed in the main part 
of the paper. 

We assume H — 70kms _1 Mpc -1 , Q m = 0.3, 
— 0.7. This only affects the Lya luminosi- 
ties and the derived star-formation rates. The re- 
ported magnitudes are on the Vega system, with 
the exception of Fig. GO 

The reduced data from this work will be avail- 
able from ESC@ and from the TOUGH websitdl 



http: / /archive, eso.org/ 



http: / /www.dark-cosmology.dk/TOUGH 
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2. Target selection, observations and data 
reduction 

2.1. The TOUGH sample 

This work is based on a parent sample named 
The Optically Unbiased GRB Host (TOUGH) 
sample. This sample of 69 Swift GRBs has sev- 
eral important features: (i) The selection criteria 
(see below) are designed to provide an optically 
unbiased (X-ray selected) sample of long-duration 
GRBs; (ii) The selection criteria are also designed 
to increase the prospects of prompt follow-up ob- 
servations being successful; (iii) The sample has 
been the focus of an extensi ve prompt follow-u p 
campaign by our group (e.g. iFvnbo et al. 2009); 
(iv) The sample has been the focus of an ex- 
tensive late-time follow-up campaign targeting 
the hos t galaxies, as reported in this series of 
papers dHiorth et al.| 120121: iMalesani et all l2012t 



Jako bsson et al. 2012t this p aper: Kriihler et al 
2012HMichalowski et al.ll2012l ). 



The sample selection cri teria and the ir ratio - 
nale are given in detail in Hiorth et al.l (2012). 
They can be summarised as follows: (1) The 
burst should trigger the 7-ray imager BAT on- 
board Swift; (2) Only long-duration bursts are 
considered; (3) An X-ray afterglow should be de- 
tected and the Swift XRT X-ray position should 
be made available within 12 hours from the trig- 
ger; (4) Milky Way extinction Ay < 0.5 mag; 
(5) Sun distance at the time of the GRB detection 
> 55°; (6) No nearby bright stars (would compli- 
cate host galaxy observations); (7) Only bursts in 
the period 2005 March 1 to 2007 August 10 are 
considered; (8) Declination in the range —70° to 
+27° (suitable for VLT observations); (9) The lo- 
calization of the burst from the X-ray afterglow 
should be better than 2.0" (90% error radius 

Furthermore, observations targeting the host as 
part of the TOUGH large program should be car- 
ried out at least 50 days after the GRB. 



6 This includes using the revised UVOT-enhanced Swift- 
XRT positions tevansl201ll;lEvans fc Osbornel20~Llj) . which 
has had the effect of in creasing the sample size from 68 (e.g. 
I Jakobsson et al.ll201lft to 69, with the additional burst be- 
ing GRB 060923B. 



2.2. Target selection 

The GRBs for the Lya host galaxy spec- 
troscopy studied here were selected from the 
TOUGH sample by applying the criterion that 
the (spectroscopic) redshift should be known and 
in the range 1.8 to 4.5. The lower limit comes 
from the atmospheric cut-off and the sensitivity 
curve of the used CCDs, while the upper limit 
comes from fringing in the CCDs used in some of 
the observing runs. 

At the time of the target selection for the last 
run of the observing campaign for the Lya spec- 
troscopy f ^2~3l . the redshift status of the TOUGH 
sample was as follows: (a) 20 bursts met the z 
= 1.8-4.5 criterion, and these were the ones ob- 
served, as listed in Table [TJ (b) 21 bursts had z 
outside the range 1.8-4.5. (c) 28 bursts did not 
have a secure, spectroscopic redshift determina- 
tion. Note that group (b) included 5 redshifts ob- 
tained as part of the T OUGH redshift campaign 



( Jakobsson et al. 20121) which were available be- 



fore the Lya observing campaign ended. 

For reference, the redshifts were subsequently 
revised for some bursts, and redshifts became 
available for other bursts. As of February 2012, 
the split of the TOUGH sample into the three 
groups based on redshift would be: (a') 27 hosts 
have z = 1.8-4.5, of which 20 hosts are those ob- 
served with FORS1 in the Lya campaign and pre- 
sented in this paper (Table [T|), while 7 hosts do 
not have such FORS1 spectroscopy^], (b') 22 hosts 
have z outside 1.8-4.5. (c') 20 hosts do not have 
a (secure, spectroscopic) redshift. 

The redshifts of two of the hosts included 
in our Lya campaign warrant special mention. 
GRB 060604 was originally included because it 
had a tenta tive afterglow redshi f t of z — 2.68 
proposed bv ICastro-Tirado et al. ( 20061 ). A sub- 
se quent re-reduc t ion an d analysis of the same data 
bv lFvnbo et~all {2009) did not confirm that red- 
shift. Instead, an upper limit of z < 3 was de- 
rived, and a possible redshift of z = 2.124 was 
suggested, based on a single absorption line in- 
terpreted as Al II. We recently obtained a n X- 
shooter host spectrum (jKruhler et al.l 120121 ) that 
gives z fa 2.1359 from H/3, [O III] and Ha. This 



7 These 7 additional bursts have redshifts from recent 
X-shooter host spectroscopy jKruhler et al.ir2012i . see also 
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Table 1 

grb sample and log of fors1 lyct host galaxy observations 



Name 


z 


Ref. 


-Rhost 

(mag) 


Grism+filtcr 


CCD 


ry-itotal 
exp 

(hr) 


Seeing 
(arcsec) 


A v 
(mag) 


GRB 050315 


1.9500 


(1) 


24.4 


600B 


new 


1.5 


0.84 


0.159 


GRB 050401 


2.8983 


(2); (8) 


26.1 


600B 


old 


2.1 


0.76 


0.216 


GRB 050730 


3.96855 


(3) 


> 27.2 


600V+GG435 


old 


1.8 


0.77 


0.168 


GRB 050820A 


2.61469 


(4) 


24.8 


600B 


old 


2.6 


0.86 


0.147 


GRB 050908 


3.3467 


(2) 


27.7 


600B 


new 


2.2 


< 1.1 


0.083 


GRB 050922C 


2.1992 


(5) 


> 26.3 


600B 


old 


2.2 


< 1.3 


0.332 


GRB 060115 


3.5328 


(2) 


27.1 


600V+GG435 


old 


2.1 


1.23 


0.447 


GRB 060526 


3.2213 


(2); (9) 


> 27.0 


600B 


new 


2.2 


0.93 


0.221 


GRB 060604 


2.136 


(6) 


25.5 


600B 


new 


1.7 


< 1.1 


0.142 


GRB 060605 


3.773 


(7) 


> 26.5 


600V+GG435 


new 


1.4 


< 1.5 


0.164 


GRB 060607A 


3.0749 


(2); (10) 


> 27.9 


300V 


new 


2.2 


< 1.4 


0.096 


GRB 060707 


3.4240 


(2) 


24.9 


600V+GG435 


new 


1.4 


1.02 


0.071 


GRB 060714 


2.7108 


(2) 


26.4 


300V 


new 


1.5 


0.99 


0.261 


GRB 060908 


1.8836 


(11) 


25.5 


600B 


new 


1.5 


< 0.9 


0.099 


GRB 061110B 


3.4344 


(2) 


26.0 


600B 


new 


2.2 


< 1.0 


0.127 


GRB 070110 


2.3521 


(2) 


25.0 


600B 


new 


1.5 


< 1.3 


0.048 


GRB 070506 


2.3090 


(2) 


26.1 


600B 


new 


1.5 


< 1.2 


0.130 


GRB 070611 


2.0394 


(2) 


> 27.0 


600B 


new 


3.0 


0.88 


0.042 


GRB 070721B 


3.6298 


(2) 


27.5 


300V 


new 


3.8 


< 0.8 


0.105 


GRB 070802 


2.4541 


(2); (12) 


25.1 


600B 


new 


1.5 


0.78 


0.090 



Note. — -Rhost is the _ff-ban d total magnitude (or 3a upper limit) of the host galaxy (before correcting 
for Galactic extinction) from iMalesani et alJ (|2012l ). CCD indicates which FORS1 CCD was used (cf. 
i )2.3p . The seeing was measured using a Gaussian fit to stars that happened to be in the slit in the 
combined spectrum. If no stars were available an upper limit on the seeing was set as the s ize of the 
smallest galaxy in the slit. Ay is the Galactic extinction in the V-band from ISchlegel et ail (| 19981 ). as 
obtained from NED. The corresponding reddening is E(B — V) = Av/3.315, and the Galactic extinction 
in the R-band is A R = 2.673 E(B - V). 



References. 



(1) iBereer et al] d2005l); (2) iFvnbo et all (|2009|); ( 3) IChen et all (|2005l ); 
(4) iProchaska et al.1 <|2007h: ( 5 )|Piranomonte et all (120081 ) ; (6)lFvnbo et alJ (12009]) , but c orrecte d for error 

Thone et all (|20ld ): (10) iFox et al.l 
redshift from this work (see text); 



see t ext); (7 ) iFerrero et al.l ll2009h : (8) 



20081); (11) iFynbo et al.l (bo09h . but 



(12) lElfasdottir et al.l (|2009h 



Watson et al l (|2006h : (9) 
prompted by the Lya 
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prompted the discovery of a wavelength calibra- 
tion error in the afterglow spectrum; the revised 
absorption-line redshift is z w 2.1361. We will 
adopt the value z = 2.136. GRB 060908 was 
originally included becau se it had an aft erglow 
redshift of z = 2.43 from iRol et all (j2006l ). Our 
spectroscopy ( §2.3[) gave a Lya host emission red- 
shift of z — 1.887. This prompted a re-analysis 
of the afterglow spectrum which did not find evi- 
dence for z = 2.43 but which did find an afterglow 
redshift of z — 1.8836 reported by iFvnbo et al 



matching our Lya host redshift 
For the target selection for the TOUGH Lya 
campaign there was no requirement that the 
hosts should be detected in the deep i?-band 
i maging from the T OUGH imaging campaign 
(jMalesani et al.l I2012T) or elsewhere. The statis- 
tics for the i?-band imaging of the 20 observed 
systems with a secure redshift in the range 1.8- 
4.5 (Table [1} are: 14 hosts are detected in the 
i?-band (with i?-band magnitudes in the range 
24.4 to 27.7) and 6 hosts are not detected down 
to a typical 3u limit of R = 27. 

All observed bursts have a detected optical af- 
terglow. This was not required, but is a con- 
sequence of the requirement of a known redshift 
before the end of the observing campaign. In 
all cases the redshift from the optical afterglow 
comes from interstellar absorption lines, both low- 
ionization metal lines (such as O I, Si II, and C II) 
and high-ionization lines (such as C IV and Si IV), 
providing a good estimate of the systemic redshift 
of the host galaxy. This is relevant for the inter- 
pretation of the velocity offset of the Lya emission 
line with respect to the afterglow redshift ( §3.21) . 

2.3. Observations 

Spectroscopic observations w ere completed us- 
ing th e FORS1 spectrograph fcf. lAppenzeller et al. 
1998) on the Very Large Telescope (VLT) over the 
period May 2006 - May 2008 in service mode. 

The detector system of FORS1 was changed 
in April 2007, i.e. during the observing campaign 
(cf. Table [T|). The old system consisted of a 
single Tektronix CCD, providing a pixel scale of 
0.20" px _1 . The new system consists of two blue- 
optimized E2V CCDs, providing a pixel scale of 
0.25" px" 1 when read out using the default 2x2 
binning as we did. The two CCDs are mounted 



so that the small gap between them is in the 
spatial direction; the gap has no practical conse- 
quences for our program. Compared to the old 
detector system, the new detector system provides 
a larger recorded wavelength range, a higher ef- 
ficiency below 6000 A, and suffers from fringing 
above 6500 A. 

All targets were observed using a 1.3" wide 
longslit. For most of the observing campaign, 
grisms 600B and 600V were used depending on the 
redshift of the target (see Table [1]) . Towards the 
end of the observing campaign the lower resolu- 
tion but higher throughput 300V grism was used 
instead of 600B for some targets. The achieved 
wavelength range and spectral resolution for the 
different grisms and detector systems are listed in 
Tabled 

The targets, which were generally too faint to 
be seen in an acquisition image, were put in the slit 
using one of two methods, both involving a nearby 
reference star. Either the position angle of the slit 
was set so that the slit would go through the refer- 
ence star and the target, or the reference star was 
put in the slit, after which an offset was applied to 
the telescope to put the target in the slit. The re- 
quired position angle or offset was computed based 
on the i?-band detection of the host or, for the 
hosts that were undetected or only marginally de- 
tected in our i?-band host imaging, based on the 
position of the afterglow. 

Each target was observed for a total exposure 
time of 1.4-3.8 hr (see Table [J), split into 4-8 indi- 
vidual exposures. The individual exposures were 
dithered along the slit. GRB 060714 was observed 
twice, since the first observation was obtained in 
poor transparency conditions. In the first obser- 
vation the galaxy continuum was not detected, 
whereas in the second observation it was. The first 
observation was done using grism 600B, whereas 
the second observation was done using the more 
efficient grism 300V and with a slightly larger ex- 
posure time. In the analysis we will only use the 
data from the second observation. 

Additionally, three bursts not in the TOUGH 
sample were observed. These are discussed sepa- 
rately in the Appendix. 
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Table 2 

Wavelength range and spectral resolution 



Grism+filtcr 


Slit width 


Wavelength range 


FWHM 


R 


c/R 


"Dispersion 






old CCD 


new CCDs 








old CCD 


new CCDs 




(arcsec) 


(A) 


(A) 


(A) 




(kms^ 1 ) 


(Apx- 1 ) 


(Apx- 1 ) 


600B 


1.3 


3300-5680 


3030-6020 


6.5 


690 


430 


1.17 


1.47 


600V+GG435 


1.3 


4510-6850 


4250-7190 


6.2 


920 


330 


1.15 


1.44 


300V 


1.3 




3000-8880 a 


14.1 


420 


710 




3.15 


600R+GG435 b 


1.3 


5090-7220 




6.5 


950 


320 


1.05 





a When grism 300V as here is used without an order sorter filter, a second order spectrum may be present for 
A > 6600 A; this has no consequences for our program. 

b Grism 600R is not relevant for the main part of this paper, only for the Appendix. 

Note. — The FWHM values were measured from the [O I]A5578A skyline in the combined science frames. The 
resolving power R — A/FWHM values were computed at the central wavelength of the available wavelength range. 
Note that the listed values of FWHM and R only correspond to the spectral resolution of the obtained host galaxy 
spectrum if the observed spatial profile of the galaxy (i.e. the intrinsic profile convolved with the seeing) was flat 
over the slit. For a peaked spatial profile the obtained spectrum will have a better spectral resolution (i.e. smaller 
FWHM and larger 7?). 



2.4. Data reduction 

Data reduction was performed mainly using 
irapH. The individual frames were bias sub- 
tracted and flat fielded. Cosmic ray events w ere 
removed using LACosmic (jvan Dokkumll200lh . A 
2D wavelength calibration was established for each 
grism and observing night (based on arc frames 
obtained the following day) and applied to the cor- 
responding science frames. Vacuum wavelengths 
were used. The wavelength calibration was ver- 
ified using the few strong skylines available, in 
particular [O I]A5578.89. The shifts in pixels in 
the spatial direction between the individual sci- 
ence frames were computed in two ways: simply 
using the requested shifts stated in the so-called 
observation blocks (OBs), and using other objects 
(preferably stars) on the slits. The two deriva- 
tions of shifts in all cases agreed to within 1 pixel 
(« 0.2"). 

The individual wavelength calibrated science 
frames were shifted in the spatial direction and 



'IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



combined (averaged). Usually no weights were 
used, but for GRB 050820A and GRB 060604, the 
signal in other objects in the spectrum varied in a 
way indicating a variable transparency, and here 
weights were used. For GRB 050820A the weights 
were set to the flux of the reference star which was 
in the slit; the weights were in the range 0.25- 
1.00 when normalized to the largest value. For 
GRB 060604 no reference star was centered in the 
slit. Based on other objects in the slit, weights in 
the range 0.67-1.00 were assigned. 

To represent the uncertainties for each pixel in 
the 2D science spectra, we calculated 2D sigma 
spectra based on photon noise and readout noise. 
The error (in ADU) in the given pixel is given by 



(1) 



where / are the counts in ADU in the given pixel 
in the 2D science spectrum before sky subtraction, 
g is the gain (conversion factor) in e~/ADU for a 
single image, n is the number of single images that 
were averaged in the combination, and RON is the 
read-out noise (in ADU) for a single image. 

The spectra and sigma spectra were flux cali- 
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brated based on sensitivity functions derived from 
30 st andard star o b servations a nd re feren c e dat a 
from iHamuy et all (| 19921 Il994h and [Ok] (|l990h . 
The spectra and sigma spectra were corrected for 
atmospheric extinc tion. The extinction curve for 
LaSilla was used (|Tusdll977l ; ISchwarz fc Melnick 
1993), since no extinction curve was available 
for Paranal at the time of the reduction. A 
comparison between the La Silla curve and the 
Paran al FORS1 broad-band extinction coeffi- 
cients IPatat 20031) shows a very good agreement 
(|Milvang- Jensen et al.ll2008l ). 

The spectra and sigma spe ctra were finally cor- 
rected for Galac tic extinction ( Cardelli et al.fll989t 
Q'DonnelllH993) using Ry = 3.1. with E{B - V) 



taken from lSchlegel et al.l ([19981 ) . Table [T] lists the 



used values. 

It should be noted that the requested observing 
conditions were to have a transparency that was 
better than or equal to thin cirrus. This means 
that the spectra were not necessarily obtained un- 
der photometric conditions, and thus that some of 
the derived fluxes may be affected by thin cirrus. 
Any such effect is mitigated by the rescaling of the 
spectra based on the photometry ( 33.1[) . 

2.5. Subtraction of neighboring objects in 
the spectra 

The spectra of two of the hosts were sub- 
stantially contaminated by a neighboring object, 
which we fitted and subtracted in the 2D spectra. 
For GRB070721B [z = 3.6) a foreground galaxy 
(z — 3.1) 1" away was fitted as a Gaussian in the 
spatial direction and a polynomial in the wave- 
length direction. For GRB 070802 the wings of a 
bright star in the slit 18" away was fitted as a poly- 
nomial in both directions. The fits were perf ormed 
using MPFIT (lMarkwardtll2009l lMorejll978h . Fig- 
ure[T]shows the spatial profiles before and after the 
subtraction of the neighboring objects. The sub- 
traction makes the derived continuum flux den- 
sities of the hosts be 3.2 times (GRB070721B) 
and 2.6 times (GRB 070802) smaller and much 
more in line with what the photometry predicts 
(cf. below). The EWs of Lya (or upper limits 
thereof) become larger by the same factors. The 
Lya fluxes, which are continuum-subtracted (cf. 
§3.1|) . are practically unaffected. 



3. Results 

3.1. Lya detections and upper limits 

We first describe the measurement of contin- 
uum flux densities in the spectra. We then com- 
pare these with the photometry and derive an av- 
erage correction for slit loss and extraction aper- 
ture loss. We finally describe the measurement of 
Lya fluxes and EWs from the spectra. 

To measure the continuum flux density in the 
spectra, apertures on the blue and red side of Lya 
were defined as follows. In rest-frame wavelength 
the apertures were 175 A wide and located such 
that a guard interval of ±6 A (corresponding to 
±1500 km s -1 ) centered on Lya placed at the af- 
terglow redshift was excluded. In the spatial direc- 
tion the continuum apertures coincided with the 
Lya apertures (see below). The width of 175 A 
was the maximum value that was covered by all 
spectra, and this value also allowed an accurate 
measurement of the continuum flux density. Un- 
certainties on all measured fluxes and flux densi- 
ties were calculated by propagating the individual 
uncertainties from the 2D sigma spectra fi j2.4[) . 

In Fig.[2]Ja) we compare the spectroscopy-based 
blue and red continuum flux densities. Of the 20 
hosts in the sample, the red continuum is detected 
for 14 hosts, while the blue continuum is only de- 
tected for 9 hosts, all at > 3<7 confidence; note 
that two of the blue non-detections are outside the 
plotted range in Fig.^a). One host (GRB 050908, 
z = 3.3, the magenta lower limit) is just above 
the 3er detection limit the blue (3.5<r) but just be- 
low it in the red (2.2c), which is plausible given 
the redshift and the sensitivity curve of the used 
grism (600B). The seemingly discrepant point at 
z = 2.3 (GRB 070506) is due to a time- variable 
pattern in the bias which we were unable to fully 
remove or quantify in terms of the error bars. This 
only has a noteworthy effect in the far blue where 
the sensitivity is very low and where the flux cali- 
bration therefore corresponds to a large amplifica- 
tion. This issue has no effect on the reported Lya 
line properties, as we have adopted the continuum 
flux densities measured in the red window in the 
spectra as those used to subtract the (small) con- 
tinuum contribution from the measured flux in the 
Lya aperture and to calculate the EWs. 

In order to investigate the absolute flux scale 
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Fig. 1— Spatial profiles of the GRB070721B and GRB 070802 spectra before (1st panel) and after (2nd 
panel) a neighboring object was fitted and subtracted in the 2D spectrum. The red thick curve shows the 
continuum on the red side of Lya, and the green curve shows Lya. The subtraction of the neighboring 
object is evident from the red thick curve. The rest of the features of this plot are explained in the caption 
of Fig. E 



T3 
<D 

O 
CD 
Q_ 



CD 

o 

CD 
Q_ 




20 40 60 80 
F x (phot,red) 



2.0 2.5 3.0 3.5 



Fig. 2. — Comparison of continuum flux densities. Panel (a) compares the continuum flux densities measured 
in the spectra on the blue and red side of Lya. Panels (b) and (c) compare the continuum flux densities 
measured in the spectra with those from the photometry (filled points: /3 = —1.5, open points [panel c only]: 
j3 = —1.0, cf. Eq. HJ, both corresponding to the red side of Lya. The units of the flux densities in panel (b) 
are 10~ 20 erg cm -2 s _1 A . The blue and red windows are centered at 1122 A and 1309 A and have widths 
of 175 A, all rest-frame. The two systems where a neighboring object was fitted and subtracted are marked: 
GRB070721B (box) and GRB 070802 (pentagon). The discrepant point in panel (a) at z = 2.3 is due to 
the blue window having very little signal and being dominated by a systematic error from a time-variable 
pattern in the bias, which has no effect on the derived Lya properties presented in this paper as they are 
based on the continuum measured in the red window in the spectra (see text). The symbol colors in the 
panels simply reflect the symbol type: red = detections, blue = upper limits, magenta = lower limits, green 
= upper limits in both the x and y direction. The spectroscopic fluxes in this figure have not been multiplied 
by a factor of 2.0 to correct for slit losses and aperture losses. 
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of the fluxes and flux densities extracted from the 
spectra, we use our F ORS2 i?-band host imaging 



(jMalesani et al.ll2012l ) . The imaging was obtained 
largely under photometric conditions, and the 



photometry was calibrated using lLandolti (|l992n . 
The derived magnitudes are total magnitudes, ob- 
tained either from using a large aperture, or from 
using a smaller aperture combined with an aper- 
ture correction. The correction was computed by 
analyzing the brighter host galaxies, and adopting 
as uncertainty the observed scatter. 

Starting from the total i?-band host magnitude 
fihost, we correct for Galactic extinction (An) to 
obtain a flux density at the observed-frame wave- 
length of the i?-band of 



Fxm = f. 



Vega 1f) 



-0.4(i? bo8t -A R ) 



where 



/^!% = 2.15xl(r 9 ergcm- 2 s- 1 . 



(2) 



(3) 



is the conversion factor for Cousins R from 



Fukugita et al.l (|1995T ); practic ally the same fac- 
tor w ould be obtained from iBlanton fc Roweisl 
(l2007lM We then extrapolate this flux den- 
sity from the effective wavelen gth of the i?-band 
(6410 A, iFukugita et"ai"1 Il995l) to the observed- 
frame center wavelength of our red window, 
(1 + z) 1309 A, by assuming an F\ oc A' 3 spec- 
trum, giving 



F x (phot, red) = F x ,r X 



(1 + z) 1309A 
6410 A 



, (4) 



where is the 
We w ill use 



rest-frame UV spectral slope. 
1 .5 as a representative value 



(e.g.. IShaplev et al.l 120031 iDouglas et ail l20ld 



201 3 k Figure Etc) shows the 



Finkelstein et al.l 
effect of as function of redshift, since in addition 
to the adopted value of = —1.5 (filled points), 
the case of = —1.0 is illustrated (open points). 

The spectroscopic flux densities show a good 
correlation with those from the photometry. This 
is shown in Fig. [5Jb) , where the two are plotted 
against each other, and in Fig. Etc) where their ra- 
tio is plotted against redshift. 12 hosts have both 
a continuum detection (at 3er) in the spectra on 



9 For Bessel R, IBlanton fc Roweij find m AB - 

m Vega = 0-21 mag and A e ff = 6442 A, which corresponds to 



the red side of Lya and a detection in the i?-band 
host imaging (at 2a). These hosts are shown as red 
points in Fig. Etc) . 2 hosts have a continuum de- 
tection in the spectra but not in the imaging; these 
are GRB 060526 and GRB 060605 and are shown 
as magenta lower limits in Fig. Etc). 2 hosts con- 
versely do not have a continuum detection in the 
spectra but have a detection in the imaging; these 
are GRB 050908 and GRB 060115 and are shown 
as blue upper limits in Fig.Et c )- 

The ratio of spectroscopic to photometric flux 
density shown in Fig. Etc) has a median value 
of 0.50 for the photometry extrapolated using 
= —1.5. We attribute the fact that this me- 
dian ratio is lower than 1 to slit losses (i.e. flux 
falling outside the slit) and extraction aperture 
losses (i.e. flux falling outside the used extraction 
apertures, as defined below), as well as possibly 
a small amount of thin cirrus affecting the obser- 
vations ( H2.4p . We use this result to derive an 
approximate global scaling factor of 1/0.50 = 2.0 
that we apply to all the flux calibrated science and 
sigma spectra ( ^2.4j) . This has the effect of mak- 
ing the derived Lya fluxes (or upper limits) and 
their uncertainties a factor of 2.0 lar ger , while the 
EWs are unaffected by this procedure! 10 !. This fac- 
tor is used throughout the paper, except in Figs. [2] 
and [8] 

The Lya fluxes were measured from the 2D 
spectra using a rectangular aperture defined in 
terms of v, the rest- frame velocity with respect 
to the afterglow redshift (cf. Table [J), and s, the 
spatial offset along the slit with respect to the af- 
terglow position. If Lya was detected, the aper- 
ture was centered on the line and the width of the 
aperture was adjusted accordingly (as illustrated 
in the last 7 panels of Figs. I3H6]). If Lya was not 
detected, an aperture of width 900 kms x 1.2" 
was used, centered at 300 kms -1 in v (a value typ- 
ical for the detected lines, cf. ^ 33 . 2j) and at s = 0.0". 
The aperture was defined in terms of integer pixels 
for simplicity. The centers and widths of the aper- 
tures are listed in Table [3l The width in v can be 
compared to the spectral resolution expressed as 
a rest-frame velocity, c/R, which is < 700 km s -1 , 
cf. Table El 



a conversion factor (Eq.(3j of 2.16x10 9 erg cm 



-1 



-1 



10 In reality the Lyc* EWs may in some cases be affected by 
slit loss, namely if the Lyo and continuum emission differ 
in term of spatial distribution. We do not have the required 
Lya narrow-band imaging to investigate this issue. 
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Fig. 3. — Lya fluxes, luminosities and (rest-frame) EWs for the 20 hosts in the sample. The dotted lines 
in panel (c) represent hosts for which no limit could be placed on the EW, due to detecting neither Lya 
nor the continuum in the spectra. The two systems where a neighboring object was fitted and subtracted 
are marked: GRB070721B (box) and GRB 070802 (pentagon). The symbol colors simply reflect the symbol 
type: red = detections, blue = upper limits. 



The fluxes within the Lya apertures in the 2D 
spectra were integrated by summing the flux den- 
sities (in erg s _1 cm~ 2 A ) and multiplying by the 
spectral pixel size (in A). Uncertainties were cal- 
culated by propagating the individual uncertain- 
ties from the 2D sigma spectra ( j)2.4p . The con- 
tinuum contribution was subtracted (if positive) 
and the errors propagated, to give the continuum- 
subtracted Lya fluxes listed in Table If the Lya 
emission was not detected at 3<r (with the a be- 
ing the uncertainty on the continuum-subtracted 
Lya flux), we give the 3<r upper limit in the ta- 
ble. Lya emission was detected in 7 hosts out 
of 20: GRB 050315, GRB 060605, GRB 060707, 
GRB 060908, GRB 070110, GRB 070506, and 
GRB 070721B. Of these, only one was known to be 
a Lya emitter fro m the literature (GRB 070110, 



Fvnbo et al. 20091) . Lya emission was not de- 



tected at 3cr for GRB 060714, only at 2.5a, but 
Lya emission was convincingly detected in the af- 
terglow spectrum, cf. £13.31 below. Note that the 
20 hosts in Figs. HHHJare sorted by Lya detection 
significance, as listed on the panels. 

Rest-frame EWs of Lya were calculated from 
the Lya fluxes and from the continuum flux den- 
sities measured in the spectra in a 175 A [rest- 
frame] wide window on the red side of Lya cen- 
tered at 1309 A (cf. above), and the errors prop- 
agated. Throughout the remainder of this paper 
it is implicit that the EWs are rest-frame values. 



The results are summarized in Fig. [3] which shows 
the Lya fluxes, luminosities and EWs (or upper 
limits in case of non-detections) versus redshift. 

The spectra of the 20 systems are illustrated in 
three figures: (i) Smoothed 2D spectra centered 
on the Lya lines are shown in Fig. 2] with the used 
aperture indicated; (ii) ID spectra (not smoothed) 
are shown in Fig. [5] and (iii) Spatial profiles (also 
not smoothed) are shown in Fig. [6] 

3.2. Velocity offset between Lya emission 
and low-ionization absorption lines in 
the afterglow spectra 

For the 7 hosts with a > 3c detection of Lya 
emission, the centroid in v within the Lya aper- 
ture was measured (see Table [3]). Since the after- 
glow redshift defines the zero point of v, this ve- 
locity measurement is the velocity offset between 
the Lya emission centroid and low-ionization in- 
terstellar absorption in the GRB afterglow spec- 
trum. A histogram of this offset for the 7 hosts is 
given in Fig. [71 The range is 200-700 km s" 1 , con- 
sistent with the few values for GRB hosts reported 
in the literature (cf. Table EJ. The histogram 
resembles the distribution of velocity offsets be- 
tween Lya emission and low-ionization interstel- 



lar absorption in L BGs (lAdelberger et al. 1200 



Shapley et al. 20031 : see also Pettini et al.l 1200 



The distribution is also in agreement with the 
velocity offsets for two Lya-selected galaxies re- 
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Table 3 

Lycc measurements from the spectra 



Name 




Lva anerture 




F(hya) 


L(Lya) 


F A (cont.) 


EW(Lya) 


u(Lya) 


c(v) 


c(s) 


w(v) 


w(s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


GRB 050315 


290 


0.08 


1102 


1.50 


23.4 ± 6.8 


0.64 ± 0.19 


86.0 ± 6.0 


9.2 ± 2.8 


283 ± 62 


GRB 050401 


326 


—0.10 


888 


1.20 


< 11.2 


< 0.80 


18.0 ± 3.6 


< 16.3 




GRB 050730 


302 


-0.04 


856 


1.20 


< 8.7 


< 1.32 


< 8.0 






GRB 050820A 


309 


0.06 


878 


1.20 


< 10.1 


< 0.56 


35.1 ±3.4 


< 8.1 




GRB 050908 


281 


-0.08 


914 


1.25 


< 6.4 


< 0.64 


< 7.0 






GRB 050922C 


342 


-0.10 


902 


1.20 


< 19.3 


< 0.70 


< 18.6 






GRB 060115 


303 


-0.10 


938 


1.20 


< 11.8 


< 1.35 


< 10.5 






GRB 060526 


294 


0.00 


941 


1.25 


< 7.9 


< 0.73 


9.8 ±2.8 


< 19.5 




GRB 060604 


307 


0.00 


806 


1.25 


< 9.0 


< 0.31 


23.9 ±3.3 


< 12.1 




GRB 060605 


564 


-0.25 


594 


1.25 


17.0 ± 2.7 


2.28 ± 0.36 


10.6 ± 2.8 


33.7 ± 10.5 


620 ± 26 


GRB 060607A 


337 


-0.05 


953 


1.25 


< 7.3 


< 0.60 


< 7.3 






GRB 060707 


788 


-1.10 


961 


2.00 


16.5 ±3.1 


1.75 ±0.33 


33.4 ±3.1 


11.2 ± 2.3 


742 ± 38 


GRB 060714 


373 


-0.05 


837 


1.25 


< 8.1 


< 0.49 


9.9 ±2.6 


< 26.3 




GRB 060908 


372 


1.00 


1002 


1.75 


77.8 ±9.5 


1.94 ±0.24 


66.8 ± 7.5 


40.4 ±6.7 


347 ±31 


GRB 061110B 


294 


0.12 


896 


1.25 


< 6.6 


< 0.71 


14.4 ±2.6 


< 10.7 




GRB 070110 


345 


-0.22 


1078 


1.75 


40.2 ±4.0 


1.73 ±0.17 


37.8 ±3.5 


31.8 ±4.3 


358 ± 26 


GRB 070506 


379 


-0.12 


983 


1.25 


13.9 ±3.5 


0.57 ±0.14 


13.0 ±3.5 


32.3 ± 11.8 


360 ± 62 


GRB 070611 


319 


-0.05 


951 


1.25 


< 9.7 


< 0.29 


< 8.7 






GRB 07072 IB 


271 


0.08 


839 


1.50 


11.2 ± 1.6 


1.37 ±0.19 


7.4 ± 1.5" 


32.5 ± 8. a 


212 ±31 


GRB 070802 


321 


0.12 


941 


1.25 


< 8.9 


< 0.43 


32.2 ± 3.1 a 


< 8.1 a 





a The listed error reflects the random error only. A systematic error due to the fitting and subtraction of a 
neighboring object is likely present. 

Note. — Columns (2)-(5) define the aperture in the 2D spectrum within which the Lya flux was measured. 
Specifically, the columns give the center (c) and width (w) of the aperture in terms of v (in kms -1 ), the rest- 
frame velocity with respect to the afterglow redshift (cf. Table [TJ, and s (in arcsec), the spatial offset along the slit 
with respect to the afterglow position. The subsequent columns are: (6) F(Lya), the Lya emission line flux, in 
units of 10 -18 erg cm -2 s -1 ; (7) L(Lya), the Lya luminosity, in units of 10 42 ergs -1 ; (8) i"\(cont.), the continuum 
flux density measured in a 175 A rest-frame wide aperture (centered at 1309 A) on the red side of Lya, in units 
of 10 -20 ergcm -2 s -1 A ; (9) EW(Lya), the rest-frame Lya equivalent width, in A; (10) u(Lya), the rest-frame 
velocity of the Lya emission line with respect to the afterglow redshift, i.e. the centroid of v within the aperture, 
in kms -1 . The listed error on ?;(Lya) is only based on the errors on the fluxes within the aperture, and does not 
include the error on the afterglow redshift. All upper limits are 3a. The flux calibrated science spectra and their 
uncertainties (sigma spectra) have been multiplied by a factor of 2.0 to correct for slit loss and aperture loss f £|3.1|) . 
This has made all the numbers in columns (6), (7) and (8) be larger by this factor. 
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Fig. 4. — The obtained 2D spectra, centered on where Lya is expected. The zero point for the rest-frame 
velocity is defined by the afterglow redshift listed in Table[TJ The zero point for the spatial (angular) position 
corresponds to the afterglow position. The shown sections have size 3800 kms -1 x 7.5". The spectra have 
been smoothed by a Gaussian with FWHM = 3px. The intensity cuts are the same for all the panels, in 
units of the noise in the given spectrum, allowing a visual comparison of the significance of the features in 
the different panels. The green rectangle marks the aperture within which the Lya flux and its uncertainty 
are measured; the aperture centers and widths are listed in Table [3] On the panels is stated by how many 
sigma Lya is detected; the panels are sorted by this. The red horizontal lines indicate where the continuum 
(if detected at > 3ct) is located, as defined by the spatial centroid in the aperture in which the continuum is 
measured, cf. Fig. [6] 
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ID spectra, derived by a straight sum over the spatial aperture (cf. Table [3]). The flux densities 

-l 



Fig. 5. 

F\ are in units of 10~ ls ergs^ 1 cm -2 A No smoothing has been applied. The detection significance of 
the Lya emission line is given on the panel; the panels are sorted by this. The vertical dotted lines mark 
the velocity limits of the Lya aperture, cf. Table |3l The horizontal dot-dashed (black) line simply indicates 
zero, while the horizontal dashed (red) line indicates the continuum, measured as the mean level in a 175 A 
rest-frame wide window on the red side of Lya. The plotted range in velocity corresponds to the range of 
the 2D spectra shown in Fig. 01 
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Fig. 6. — Spatial profiles, i.e. the flux densities averaged over a given wavelength range versus spatial 
coordinate in the 2D spectrum. Red thick curve: (F\) calculated over the 175 A rest-frame wide continuum 
window on the red side of Lya, in units of 10~ 20 ergs -1 cm~ 2 A . Green curve: (F\) calculated over the 
wavelength range of the Lya aperture, in units of 10~ 19 ergs -1 cm -2 A . No smoothing has been applied. 
The detection significance of the continuum in the above-mentioned red window is given on the right hand 
side of the panel in red, and the detection significance of the Lya emission line is given on the left hand side 
of the panel in green. The vertical dotted lines mark the spatial limits of the Lya aperture, cf. Table [H The 
plotted range in spatial coordinate corresponds to the range of the 2D spectra shown in Fig. 01 The panels 
are sorted by Lya detection significance. 
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ported bv lMcLinden et al.l (|201lh using the [O III] 
emission line to define the systemic velocity. It 
should be noted that the assumption that the af- 
terglow redshift provides the systemic velocity is 
only valid on average (over a sample of hosts) . For 
individual hosts the GRB sightline may probe a re- 
gion of the host that has a non-zero velocity due 
to the internal kinematics of the galaxy (e.g. rota- 
tion). This cannot be a large effect, since other- 
wise the measured velocity offsets of Lya (Fig. [7]) 
would not all have the same sign. 

It should be noted that what we measure is sim- 
ply the centroid of the Lya emission line in our 
GRB host spectra, which may not be identical to 
the peak of the line if the line is asymmetrical. 
Our spectra (Fig. [5]) do not have sufficient spec- 
tral resolution or S/N to investigate this issue. 




417 km s" 
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400 600 800 1000 
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Fig. 7. — Distribution of velocity offsets between 
the centroid of the Lya emission in the GRB host 
spectrum and low-ionization interstellar absorp- 
tion lines in the GRB afterglow spectrum. 

The origin of the offset is most likely a combi- 
nation of radiative transfer of the resonantly scat- 
tered Lya photons and a star-formation driven 
outflow from the host galaxy (for a full discus- 



sion of these effects we refer to lFvnbo et al.l l2010. 
their §4.3). It should be emphasized that the 
velocities presented in Table [3] and Fig. [7] sim- 
ply represent the centroid of the Lya line in the 
spectrum (with respect to the afterglow redshift); 



they do not directly translate into outflow veloc- 
ities. In the outflow scenario several factors af- 
fect the observed redwards shift of the Lya ve- 
locity centroid with respect to the systemic veloc- 
ity. High column densities and (to a lesser ex- 
tent) low temperatures push t he Lya peak fur - 
ther from the systemic velocity ( Harrington! 1 973h . 
The velocity shift also increases with increasing 
outflow velocities, up to ~ 10 kms , where the 
peak starts to drift back toward t he systemic ve- 
locity (e.g. Verhamme et al.ll2006l ). On the other 
hand, if dust is present, preferentially the wings of 
the line are removed , effectively reducing the shift 



( Laursen et al. 20091 ). This effect is stronger the 



more homogeneous the medium is, since dumpi- 
ness of the g as and dust facilitates the escape o f 
Lya photons (|Neufeldlll99lt lHansen fc Ohll2006|) . 



3.3. Lya emission comparison: host and 
afterglow spectra 

For three of the hosts in the sample Lya 
emission was detected directly in the afterglow 
spectra: G RB 060714 (Uakobsson et al.l l2006bh . 
GRB 070110 dFvnbo et al.l2009fl and also marginally 
in GRB 070721B (|Fvnbo etall 
compares the spectra. For 



2009). 



Figure |8] 
GRB 070110 and 
GRB 070721B the Lya flux measured from the af- 
terglow and host spectroscopy is consistent within 
the errors (at 2 sigma). However, for GRB 060714 
there is a significant difference, with the Lya flux 
from the afterglow spectrum being almost a fac- 
tor of 4 larger than in the host spectrum reported 
here. In the afterglow spectrum the Lya emission 
appears to be very extended both spatially and 
in velocity space. In addition, the position angles 
of the two observations are nearly perpendicular 
(0° east of north for the afterglow spectrum and 
116° for the host galaxy spectrum). Hence, we 
suspect the cause of the difference is a higher slit 
loss in the host galaxy spectrum. It is possible 
that other of the 12 hosts where we do not detect 
Lya emission at > 3er actually would have been 
detected as Lya emitters had we used a slit at a 
different position angle or a wider slit. 

4. Discussion and summary 

In this work we have carried out a systematic 
search for Lya photons from GRB host galax- 
ies selected from the larger well-defined TOUGH 
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Fig. 8. — Comparison of host spectra (red solid 
lines) with afterglow spectra (blue dotted+solid 
lines). The three GRBs shown are those for 
which Lya emission was detected in the after- 
glow spectrum superimposed on the damped Lya 
absorption trough. The position angle (PA) of 
the slit is given on the panels; note the large 
PA difference between host and afterglow spec- 
troscopy for GRB 060714. The spectra have 
been corrected for Galactic extinction. F\ is 
given in units of 10 -18 erg cm -2 s^ 1 A . In 
the comparison for GRB070721B it should be 
noted that only the host spectrum has had a 
neighboring galaxy subtracted (cf. §3.1|) . a pro- 
cedure that decreases the overall flux level by 
about 0.1 in the plotted unit s. The after- 
glow s pectra were taken fr om Jakobsson et al.l 
(|2006bh [GRB 060714] and iFvnbo et all (|2009h 
[GRB 070110 and GRB070721B]. None of the 
spectra in this figure have been corrected for slit 
losses or aperture losses. 



sampl e of such galaxies presented in Hiorth et al 
12012). Unlike previous studies (cf. IFvnbo et al.. 
20031 ) we find that Lya emission is not ubiquitous 
among GRB host galaxies. Of the 20 host galax- 
ies studied here we detect (at 3a) Lya emission 
from 7 of them (with the Lya [rest-frame] EW 
in the range 9-40 A), derive 3cr upper limits on 
the Lya EW for 7 of them (in the range 8-26 A), 
while we obtain no constraints on the Lya EW for 
the last 6 hosts (due to neither detecting the con- 
tinuum nor Lya emission in the spectra, both at 
3a), cf. Table El Out of the 14 hosts with either 
a Lya EW or an upper limit on the EW, 8 hosts 
have Lya EW less than 20 A (rest-frame) , which is 
the typical limit in narrow-band surveys for Lya 
emitters. For the 7 detections, the measured EWs 
in the range 9-40 A are low compared to the dis- 
tribution of EWs found for narrow-band selected 
galaxies at similar redshifts dGrqnwall et aL 2007 : 
Grove et al1l2009l; iNilsson et al.ll2009h . 

The Lya luminosities for the 7 GRB hosts 
with detected Lya emission are in the range (0.6- 
2.3) xl0 42 ergs _1 . Such fairly low Lya luminosi- 
ties are only pro bed by a few studies of Lya emit- 
ting g a laxies, e.g. Rauch et all ( 20081) ; iGrove et al.l 
(|2009l ); ICassata et al.l (|201ll ). The Lya luminosity 
can be translated into a star-formation rate (SFR), 
assuming no dust extinction, as 



SFR: 



i(Lya) 



1.1 x 10 42 ergs" 



■ Mi 



(5) 



using the r elation between SFR and L(Ha) from 
Kennicuttl (|l998l ) and the c ase B recombinatio n 
ratio L(Lya)/L(Ha) = 8.7 (jBrocklehurstl Il97ll) . 
The observed range in Lya luminosity for the 7 
detections would translate into a range in SFR of 
0.6-2.1 Mq yr _1 , but the assumption of no dust 
extinction is likely not always correct, as illus- 
trated by the trends of Lya luminosity and EW 
with afterglow spectral index discussed in H4.2I If 
dust is present then Eq. ([5]) provides a lower limit 
of the SFR. 

4.1. Comparison with LBGs 

Lyman-break selection and GRB selection are 
complementary methods to identify samples of 
galaxies at high redshift. In this section and in 
Fig.Hllwe carry out a comparison of the LBGs from 
IShaplev et al. ( 2003 ) with the GRB host galaxies 
from this work. 



16 




50 100 
EW(Lya) [A] 



150 200 



Fig. 9.— 

LBGs from 



Comparison of L ya E Ws between the 



Shaplev et al 



(l2003h and the GRB 
host galaxies from this work. Small black points 
and grey/black histogram: LBGs (N = 803). 
Other points: GRB host galaxies (N — 14), in (a) 
and (b) split into galaxies with detected Lya emis- 
sion (red filled circles and histogram) and without 
(blue open circles and histogram). 



The 20 GRB host galaxies from this work fall 
into 3 categories. For 6 galaxies we detected nei- 
ther the continuum nor Lya emission in the spec- 
tra, and these galaxies are omitted from the analy- 
sis. For 7 galaxies we detected both the continuum 
and Lya emission in the spectra, and for these we 
use the Lya EWs reported in Table[3] For another 
7 galaxies we detected the continuum but not Lya 
emission in the spectra. In Table [3] and the rest 
of the paper we have reported the 3cr upper lim- 
its on the Lya EWs. An example is GRB 061110B 
where the measured [rest-frame] EW is 5.0 ±3. 6 A, 
which we replaced by the 3tr upper limit of EW 
< 10.7 A. This procedure implicitly assumes that 
Lya can only be in emissi on. However, the LBGs 
from IShaplev et al. (2003) often show significant 
Lya absorption (negative EWs), so in order to 
make a fair comparison with that sample, we use 
the measured EWs also for the 7 GRB host galax- 
ies without detected Lya emissio 

Figure [Ha) shows apparent i?-band magni- 
tude versus Lya EW. The 803 LBGs are shown 
as small black open squares. The 7 GRB host 
galaxies with Lya emission detected at 3tr are 
shown as red filled circles, and the 7 GRB host 
galaxies without such detected Lya emission are 
shown as blue open circles. The plot shows that 
the GRB hosts from t his work typica l ly are fainter 



than the LBGs from IShaplev et all (|2003l ). This 



is also the case for the luminosities, since the red- 
shift distributions of the two samples are fairly 
similar: LBGs: (z) = 3.0, sd = 0.3; GRB host 
galaxies: (z) = 2.8, sd = 0.6 (with sd being 
the standard deviation). It is tempting to de- 
fine a faint subset of the LBG s ample that is 



better matched to our sample, but S haplev et al 



(|2003l ) conclude that the redshift incompleteness 
at fainter ma gnitudes (say fainter than 1Z ss 24.5, 
cf. Fig. 7 in IShaplev et al.1 120031 ) is likely such 
that preferentially galaxies without (strong) Lya 
emission are missing. This might argue for only 
comparing our GRB hosts with a bright LBG sub- 
sample, but then the luminosity difference would 



"These EWs are: GRB 050401: -4.4±5.4 A; GRB 050820A: 
-4.2±2.7A; GRB060526: 4.4±6.5A; GRB060604: 3.0± 
4.0 A; GRB 060714: 18.3±8.8A; GRB061110B: 5.0±3.6A; 
and GRB070802: -3.5±2.7A. 

12 The apparent magnitudes in Fig. |9fa) are on the AB sys- 
tem and have been corrected f or Galactic ext i nction . The 
filt er used for the LBGs from IShaplev et~aT1 (I2003T) is K, 
see lSteidel fc Hamilton! l|l993l 1. 
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be substantial. We will therefore simply use the 



full IShaplev et al.l (|2003l ) for comparison with our 
GRB host sample. 

Figure |^b) shows histograms of the EWs: 
grey filled histogram: LBGs; hacthed histograms: 
GRB hosts galaxies, with blue and red having the 
same meaning as in panel (a). The LBG histogram 
has been scaled down by a f actor of 50, but is oth - 
erwise identical to Fig. 8 in lShaplev et al. I (l2003h . 

Figure [9jc) shows the cumulative EW distribu- 
tions: smooth black curve: LBGs; jagged magenta 
curve: the 14 GRB host galaxies. A Kolm ogorov- 
Smirnov (K-S) test (e.g. iPress et al.lll992l) gives a 
1.7% probability that the two samples are drawn 
from the same parent distribution. In other words, 
we detect a difference at 98.3% confidence between 
the two samples. This result is driven by the lack 
of GRB host galaxies with substantial Lya absorp- 
tion (i.e. with EWs below —5 A). 

iFvnbo et al. (2003) also found a significant dif- 
ference (99.8% confidence) between the EWs of 5 
pre- Swift GRB host galaxies an d an approxima- 
tion of the IShaplev et all ( 2003 ) distribution. If 
we use our updated compilation of the EWs for 
these 5 pre -Swift hosts and carr y out a K-S test 
against the lShanlev et al.l (|2003l) sample, we get a 
similar result, namely a difference that is signifi- 
cant at 99.2% confidence. We compare the pre- 
Swift sample with the sample from this work in 



4.2. The relation between afterglow spec- 
tral index and host Lya emission, and 
comparison with pre-Swift studies 

Remarkably, substantially larger EWs were 
found in the previous, pre- Swift studies of Lya 
emission from GRB hosts despite the fact that 
these studies targ eted a much smaller sample 
(jFynbo et al.ll2003l and references therein). Our 
updated compilation of the EWs for th e 5 pre- 
Swift hosts studied bv lFvnbo et al. ( 2003 ) is given 
in the first 5 rows of Table 02 the 3 large values 
around 70-100 A are noteworthy. A K-S test com- 
paring the EWs of the pre-Swift sample (JV = 5) 
with the sample from this work (N — 14, cf. $4.1\i 
gives an 8% probability that the two samples are 
drawn from the same parent distribution. This is 
marginal evidence for a difference. This difference 
could therefore be a chance effect, but another 



plausible explanation is different biases in the two 
samples. 

The present sample is based on an underlying 
X-ray selected sampl e of 69 bursts (the TOUGH 
sample, see ^2. II and Hiorth et al. 2012?) which is 
nearly unbiased. The sample of 20 bursts fol- 
lowed up for Lya spectroscopy in this work (i.e. 
those with a known afterglow redshift in the range 
1.8-4.5) is biased since an optical afterglow was 
de facto required, and since some bursts in the 
TOUGH sample were without a determined red- 
shift at the time of the target selection for the Lya 
spectroscopy and thus could be in the targeted 
redshift range of 1.8-4.5 (indeed, 7 of these bursts 
were recently found to be at z = 1.8-4.5 from X- 
shooter host spectroscopy, see lKruhler et al" 2012 
and below, while 20 of the TOUGH bursts still do 



not have a determined redshift, cf. 



The pre- 



Swift sample of 5 bursts (jFvnbo et al.l 120031 and 
references therein) is even more biased towards rel- 
atively bright optical afterglows due to the larger 
times to localize the b urst and larger lo calization 
uncertainties (see also Kann et al. 2010h . This is 
shown in Fig. [TU^a) , where we plot the afterglow 
i?-band magnitude at 12 hr after the burst (see 
Table [5]) versus redshift for the pre-Swift sam- 
ple (open green stars) and the Swift sample from 
this work (other symbols). In panel (b) we plot 
EW(Lya) (detections, and for our sample also up- 
per limits) versus afterglow magnitude. Compar- 
ing the two samples suggests that the larger Lya 
EWs for the pre-Swift sample is related to brighter 
afterglows, which in turn could be related to galax- 
ies having less dust. On the other hand, within the 
sample from this work there is no evidence for a 
correlation between Lya EW and afterglow mag- 
nitude. 

To further examine the role of dust we turn to 
the afterglow spectral index /3ox, defined by 



Pox = 



log [i^QoptV^te)] 
log [z/x/fopt] 



(0) 



where F v is the flux density of the afterglow and 
where f op t and z^x are representative center fre- 
quencies (pivotal frequencies) of the optical and 
X-ray bands, respectively. If F v were a single 
power-law between v opt and vy^ it would have the 
form F v oc v~$ m . A low value of /3qx indicates 
suppression of the optical emission compared to 
the X-ray flux. For low-redshift events (e.g. for 
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Table 4 

Afterglow i?— band magnitudes at 12 hr after the burst — The 5 PRE-Swift Lyo. emitters 

AND THE 20 HOSTS FROM THIS WORK 



Name 


EW 


R 


Ref. 


Comment 


GRB 971214 


yes 


22.06 ± 0.06 


(1) 




GRB 000926 


yes 


18.33 ±0.10 


(2) 


Extrapolated from R = 19.33 at 20.7 hr using a — 1.69 


GRB 011211 


yes 


20.23 ± 0.04 


(3) 




GRB 021004 


yes 


18.60 ± 0.03 


(4) 




GRB 030323 


yes 


19.00 ±0.1 


(5) 




GRB 050315 


yes 


20.90 ± 0.20 


(6) 




GRB 050401 


UL 


23.00 ±0.10 


(7) 




GRB 050730 




20.20 ± 0.1 


(8) 




GRB 050820A 


UL 


18.78 ± 0.14 


(6) 




GRB 050908 




22.38 ±0.10 


(9) 


Data at 0.68, 4.25 and 22.12 hr, a = 0.78 


GRB050922C 




19.98 ±0.03 


(6) 




GRB 060115 




21.77 ±0.12 


(10) 




GRB 060526 


UL 


19.82 ±0.05 


(11) 




GRB 060604 


UL 


21.40 ±0.30 


(12) 




GRB 060605 


yes 


20.45 ±0.1 


(13) 




GRB 060607A 








Unusual light curve — unable to interpolate 


GRB 060707 


yes 


21.30 ±0.30 


(14) 


Interpolation between two datapoints 


GRB 060714 


UL 


21.10 ± 0.1 


(15) 




GRB 060908 


yes 


21.90 ± 0.1 


(16) 




GRB 061110B 


UL 


> 23.30 


(17) 


Limit assuming a — 0.5 from the VLT observation at 2.3 hr 


GRB 070110 


yes 


20.00 ±0.1 


(18) 


Converted from U-band 


GRB 070506 


yes 


> 20.14 


(19) 


Limit assuming a — 0.5 from the VLT observation at 3.6 hr 


GRB 070611 




21.50 ±0.30 


(17) 


Extrapolated from R = 21.0 at 7.7 hr adopting a = 1.0 ± 0.5 


GRB 070721B 


yes 


22.50 ±0.30 


(20) 


Extrapolated from R = 23.0 ± 0.1 at 17.6 hr adopting a = 1.0 ± 0.5 


GRB 070802 


UL 


23.50 ±0.1 


(21) 





Note. — The first 5 bursts are from the pre-Swift sample l|Fvnbo et al.l [20031 . and references therein), and the re- 
maining 20 bursts are from the Sura/f-based sample studied in this work. The EW column indicates what type of 
information is available about the Lya EW of the host (detection, upper limit or no constraint); the actual values are 
given in Tables [3] and [5] R is the afterglow /2-band magnitude at 12 hr after the burst. In the cases where we directly 
read the afterglow magnitude at 12 hr after the burst from a plot in the stated reference we have assigned a magnitude 
error of 0.1. The parameter a is the slope of the light curve: F(t) oc t~ a . 



(1) iDiercks et all Jl99St); (2) iFynbo et all |200ll); (3) [jakobsson et all J2004bl); (4 ) . I Holland et al 

20101 ): (71 Iwatson et all |2006l ); f8) |Pandev et all (|2006l ) 



Ref erences 

(|2003l ); (5) IVreeswiik et all (|200^ T~(6l iKann et al 
measurement using archival data from ESO program 275.D-5022 (PI 
arch ival data from ESQ prog ram 076.A-0392 (PI: Tagliafer ri) ; (11 
(13)| Ferrero et all |2009h : (14) "~ 
fl6) ICovino et all (^Old ); (17) 
typo: correct value is R 



Chincarin i) 



Thone et al 



de Ugarte Postigo et all ( 20061) and Jakobsson et al 
Fvnbo et all (|2009T l 



(9) Our 

(10) Our measuremen t usin g 
2oioh : (12) iTanviret all | 2006 f) : 



a ..2006ah; r i 5)lKrimm et all 

(18) iTroia et all (120071 ): figl lFvnbo et all (|2009h but corrected for 



2007); 



19.5 at 3-6 hr- 
(PI: Vreeswijk); (21) iKriihler et all |2008l 



(20) Our measurement using data from our own ESO program 079.D-0429 
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Fig. 10. — Afterglow ii-band magnitude at 12 hr after the burst, versus redshift (panel a) and Lya EW 
(panel b). Open green stars: the pie-Swift sample of 5 GRB host Lya emitters. Other symbols: the sample 
of 20 hosts from this work: red filled circles: hosts with a Lya detection, blue filled cirles: hosts without 
a Lya detection but with an upper limit on the Lya EW, blue open circles (panel a only): hosts with no 
contraints on the Lya EW. 



z < 4.5 as considered here), where the optical is 
not cut off by inter-galactic medium absorption, 
Pox is thus connect ed to dust extinctio n along the 
GRB sightline (e.g. lFvnbo et al.ll2009h . In partic- 
ular, assuming standard synchrotron theory, pox 
cannot be intrinsically smaller than 0.5, and there- 
fore bur sts with Pox < 0-5 are r eferred to as "dark 
bursts" ( Jakobsson et al.|[2004al ) . although moder- 
ate extinction can be present also in bursts with 
larger values of Pox- 

We have compiled a list of all known GRB hosts 
with Lya emission, both from this work and from 
the literature, including both pie-Swift and Swift 
bursts, see Table [5l The table lists both the after- 
glow Pox and the host Lya emission properties. 
We have corrected the literature Lya fluxes for 
Galactic extinction where needed. The table only 
contains hosts with a Lya detection. The Lya up- 
per limits from this work are in Table [3] ( with P ox 
available for all bursts from lFvnbo et al1l2009l ). 

In Fig. [Tl] we plot host EW(Lya) and host 
L(Lya) versus afterglow Pox, both for the hosts 
from this work and for the additional hosts with 
Lya detections from the literature. The plots 
show a lack of hosts with high EW(Lya) and large 



L(Lya) at the low end of the /3ox range. This indi- 
cates that dust extinction is important in reducing 
the strength of the Lya line. GRB061222A goes 
against the trend, with a very low Pox (namely 
an upper limitEl of /3 ox < 0.22) and detected Lya 
emission. This can be explained by the fact that 
Pox only probes the afterglow sightline, whereas 
the host Lya emission is a quantity that is global 
for the galaxy. There are indeed several cases 
where a dark GRB exp loded in an overall blue 
galaxy, e.g. GR B 070306 (Ijaunsen et al.l l2008l) and 
GRB100621A (|Kriihler et all 1201 ll ). It is also 
seen from Fig. [TT] that the hosts of the pie-Swift 
GRBs with high EW(Lya) and large L(Lya) are 
preferentially found at the high end of the Pox 
range. This suggests that the pie-Swift sample 
(shown as op e n gree n stars in Fig. [II]) discussed by 



Fvnbo et al.l ([20031 ) is more biased against dusty 



sightlines. 

Our finding that Lya emission is not ubiq- 



13 The upper limit comes from the afterglow only having an 
f?-ba nd upper limit. Th e afterglow has a Kg-band detec- 
tion iCenko fc Foxll2006t) . which gives /3 Q x = 0.10 (with 
"O" now signifying K s — band rather than optical//?— band), 
which is even more constraining. 
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Table 5 

Known GRB host Ly« emitters: afterglow spectral slope and host Ly« properties 



Name 


z 


T 


n 

POX 


Ret. 


hjW(Lya) 


net. 


F(Lya) 


Ret. 


L(Lya) 


v(Lya) 


Ket. 


GRB 971214 


3.42 


no 


0.64 


(1) 

V- 1 -/ 


13 


(6) 


6.2 ± 0.7 


(13) 


0.66 ± 0.07 






GRB 000926 


2.04 


no 


0.87 


(1) 


71 ± 18 


(7) 


149 ± ll b 


(7) 


4.51 ±0.33 






GRB 011211 


2.14 


no 


0.98 


(1) 


21 ± 10 


(8) 


33.6 ± 9.6 b 


(8) 


1.14 ±0.33 






GRB 021004 


2.33 


no 


0.93 


(1) 


68 ± 11.5 


(9) 


313 ± 64 b 


(14) 


13.14 ±2.67 


530 


(14) 


GRB 030323 


3.37 


no 






108 ± 38 


(10) 


12 ± 1 


(15) 


1.23 ± 0.10 


151 ±46 


(15) 


GRB 050315 


1.95 


yes 


0.63 


(2) 


9.2 ±2.8 


(11) 


23.4 ±6.8 


(11) 


0.64 ±0.19 


283 ± 62 


(11) 


GRB 060605 


3.77 


yes 


1.00 


(2) 


33.7 ± 10.5 


(11) 


17.0 ± 2.7 


(11) 


2.28 ± 0.36 


620 ± 26 


(11) 


GRB 060707 


3.42 


yes 


0.73 


(2) 


11.2 ± 2.3 


(11) 


16.5 ±3.1 


(11) 


1.75 ±0.33 


742 ± 38 


(11) 


GRB 060714 


2.71 


yes 


0.77 


(2) 






17.3 b 


(16) 


1.05 






GRB 060908 


1.88 


yes 


0.80 


(3) 


40.4 ±6.7 


(11) 


77.8 ± 9.5 


(11) 


1.94 ±0.24 


347 ± 31 


(11) 


GRB 060926 


3.21 


no 


0.87 


(4) 






62.1 ±4.9 b 


(17) 


5.65 ± 0.45 


311 


(17) 


GRB 061222A 


2.09 


no 


<0.22 


(2) 


31 


(12) 


168 b 


(18) 


5.39 






GRB 070110 


2.35 


yes 


0.77 


(2) 


31.8 ±4.3 


(11) 


40.2 ± 4.0 


(11) 


1.73 ±0.17 


358 ± 26 


(11) 


GRB 070506 


2.31 


yes 


0.93 


(2) 


32.3 ± 11.8 


(11) 


13.9 ±3.5 


(11) 


0.57 ±0.14 


360 ± 62 


(11) 


GRB 070721B 


3.63 


yes 


0.72 


(2) 


32.5 ± 8.0 a 


(11) 


11.2 ± 1.6 


(11) 


1.37 ±0.19 


212 ± 31 


(11) 


GRB 071031 


2.69 


no 


0.97 


(2) 






23.6 ± 2.7 b 


(17) 


1.41 ± 0.16 


254 


(17) 


GRB 090205 


4.65 


no 


0.98 


(5) 






23.6±4.9 b 


(19) 


5.17 ± 1.08 


180 ± 153 


(19) 



a The listed error reflects the random error only. A systematic error due to the fitting and subtraction of a neighboring 
object is likely present. 

b The published Lya flux or the provided spectrum was not corrected for Galactic extinction, but we have applied the 
correction. 



Note. — T indicates whether the host is part of the TOUGH sample ( q2.ll H iorth et al.ll2012j ) studied in this work. /3ox is 
the afterglow optical-to- X-ray spectral slope (see Eq. [6]), where optical means 7?-band (unless otherwise stated) and X-ray 
means 3 keV. EW(Lya) is the rest-frame Lya emission line EW, in A. F(Lya) is the Lya emission line flux, in units of 
10 -18 ergcm~ 2 s -1 . L(Lya) is the Lya emission line luminosity, in units of 10 42 ergs _1 . i;(Lya) is the rest-frame velocity 
centroid of the Lya emission line with respect to the afterglow redshift, in km s _1 . "Ref." gives the reference for the preceeding 
column. The bursts up to and including GRB 030323 are pie-Swift, while the remaining bursts are from Swift. All Lya fluxes 
and luminosities are corrected for Galactic extinction. Note that GRB 030429 is not included, since even though its spectrum 
showed an indication of Lya emission, it was not statistically significant (< 2a) |jakobsson et al]|2004d ) ■ 



References . - (1) IJakobsson et al.l d2004a); (2) Fvnbo et all 120091 ); (3) Our calculation, usin g _R-band data fro m 
ICovino et all l|2010Tl and X-ray data from lEvans et all (|2007i . 120091 ) (we note that the /3ox value in iFvnbo et all 120091 is 
in er ror); (4) Our calcu l ation, using archival ii-band data from ESO program 077.D-080 5 (PI: Tagl i aferri ) and X-ray data 
from [Evans et al | (|2007l . l2009h : (5) Our calculation, using 7-band photomet ry from lD'Avanzo et al.l (|2010l ) (rather than R- 
band , due to the high redsh ift) and X-ray data fromlEvans et all (|2007i . l2009l ); (6) Our measurement, using the host spectrum 
from iKulkarni et all (Il998|); (7 ) lFvnbo et all (I2002TI. with F (Lya) corrected for continuum contribution as p rescribed in that 
aper ; (8) IFvnbo et all ( 20031 '): (9) Jakobsson et all l|2005l ); (10) Our calculation, using the line flux from IVreeswiik et all 
20041 ') and deriving the continuum flux density from the F606W photometry from Wa inwright et all |2007l ) (correcting for 
Lya forest absorption and imposing a slitl oss; t hese effects almost can cel ou t); (11) This work; (12) Our measuremen t , usin g 
the host spectrum fromlPerlev et al.l (|2009l ); (13) lKulkarni et all j 19981 ); ( 14) llVtoller et all d2002l); (15)IVreeswiik et all (|2004l ); 
(16) IJakobsson et all d2006bl) ; (17) Our m easurement, using the afterglow spectrum from lFvnbo et al.l (|2009l ); (18) D. Perley, 
priv. comm.. cf. lPerlev et all l|2009D ; (19) iD'Avanzo et all l|20lb1 ). 
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Fig. 11. — Host Lya EW (panel a) and host Lya luminosity (panel b) versus afterglow spectral index /3ox- 
The GRBs in the sample of this paper are shown as red filled circles (Lya detections), blue upper limits 
(Lya non-detections), and dotted lines (unconstrained Lya EWs, panel a only). Additional GRB hosts with 
Lya detections from the literature are shown as green stars (open stars: pre-Swift, hlled stars: Swift). These 
bursts are not in the sample of this paper, except for GRB 060714: this burst has a Lya detection from 
the literature, whereas our data resulted in an upper limit, as discussed in i)3.3l This burst is therefore 
plotted twice (panel b only): as a detection and as an upper limit. Some of the literature bursts plotted 
in panel (b) are absent from panel (a) due to not having a measured EW, which in turn is due to only an 
afterglow spectrum being available. The plotted data for all the Lya detections (both from this work and 
from the lit erature) are g i ven in Table [5] The non-detections (this work only) are given in Table El with /?ox 
taken from iFvnbo et alj (|2009h . Note that GRB 030323 (listed in Table [5) is not plotted here due to not 
having X-ray observations and hence no /3ox- The two systems where a neighboring object was fitted and 
subtracted are marked: GRB070721B (box) and GRB 070802 (pentagon). Note the logarithmic axis for the 
Lya luminosity. 
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uitous among GRB host galaxies has implica- 
tions on how well GRBs trace the overall mas- 
sive star-formation activity and on the nature of 
GRB progenitors. Given that Lya photons are 
more easily destroyed by dust than other UV 
photons due to resonant scattering, it has been 
argued that GRB hosts have low dust content. 
This could be due, among other reasons, to low 
metallicity, in agree ment with the prediction of 
the collapsar model ( MacFadven fc Woosley|[ l999; 
lYoon fc Langerl 120051: IWooslev fc Hegerl 120061) . 
Our analysis of a larger sample of GRB hosts 
shows that Lya emission is less ubiquitous than 
previously found based on a much smalle r sample 
( Fvnbo et al.ll2003l; I Jakobsson et al.ll2005 ). so that 
the above argument is not valid. Whereas other 
mechanisms than dust can reduce the strength of 
the Lya line (e.g. the geometry of the interstel- 
lar medium), the trend visible in Fig. ITT1 suggests 
that the strength of the Lya line is related to the 
presence of dust. We note that the objects in 
the sample studied in this work all have a red- 
shift measured from the optical afterglow, hence 
they are biased against very dusty systems. If the 
connection between the presence of dust and the 
weakness of the Lya line holds, we expect that 
the hosts of optically-obscured (i.e. dark) GRBs 
should have even less prominent Lya emission. 

The recent work of iKruhler et all (|2012h pro- 
vides additional insight. VLT/X-shooter was used 
to target several TOUGH hosts that lacked red- 
shifts. For 7 of the TOUGH hosts the found red- 
shift was in the range z = 1.8-4.5, and these are 
thus hosts missed by the target selection for this 
work (cf. £12. 2p . The redshifts were based on de- 
tecting one or more of the following emission lines: 
[O II], H/3, [O III] and Ha. In no cases was Lya 
detected. These 7 bursts mostly have low 0ox val- 
uei 
be 



While the lack of Lya emission still has to 
quantified in te r ms of upper limits on the EWs, 



the iKriihler et al.l (|2012l ) result supports the pic- 



ture that weak or absent Lya emission is at least 
in part caused by dust. 

The Dark Cosmology Centre is funded by the 



14 The /3 x values are: GRB 050819: <0.90, GRB050915A: 
<0.44, GRB 051001: <0.56, GRB 060814: <-0.06, 
GRB 070103: <0. 48, GRB 070129: 0.62, and 
GRB070419B: 0.25 iFvnbo et al.l |2009| , and references 
therein). 



Danish National Research Foundation. BMJ and 
JPUF acknowledge support from the ERC-StG 
grant EGGS-278202. This research has made 
use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and 
Space Administration. This work made use of 
data supplied by the UK Swift Science Data Cen- 
tre at the University of Leicester. We thank 
S. George Djorgovski for providing the spectrum 
of the GRB 971214 host and Dan Perley for pro- 
viding the spectrum of the GRB061222A host, 
cf . Table [5j Peter Laursen is thanked for discus- 
sions about Lya radiative transfer. Alice E. Shap- 
ley is thanked for discussio ns and for providim 



(200 



the data on the LBGs from IShaplev et al 
in machine-readable format. Thomas Kriihler is 
thanked for discussions about the X-shooter host 
spectroscopy. We thank the anonymous referee for 
comments that helped improve the presentation. 
Facilities: VLT (FORS1). 

REFERENCES 

Adams, T. F. 1972, ApJ, 174, 439 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., 
& Pettini, M. 2003, ApJ, 584, 45 

Appenzeller, I., Fricke, K., Furtig, W., et al. 1998, 
The Messenger, 94, 1 

Berger, E., & Becker, G. 2005, GRB Coordinates 
Network, 3520, 1 

Berger, E., Kulkarni, S. R., Fox, D. B., et al. 2005, 
ApJ, 634, 501 

Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 

Brocklehurst, M. 1971, MNRAS, 153, 471 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 
1989, ApJ, 345, 245 

Cassata, P., Le Fevre, O., Garilli, B., et al. 2011, 
A&A, 525, A143+ 

Castro-Tirado, A. J., Amado, P., Negueruela, I., 
et al. 2006, GRB Coordinates Network, 5218, 1 

Cenko, S. B., & Fox, D. B. 2006, GRB Coordinates 
Network, 5978, 1 



23 



Chariot, S., & Fall, S. M. 1993, ApJ, 415, 580 

Chary, R., Bcrger, E., & Cowic, L. 2007, ApJ, 671, 
272 

Chen, H.-W., Prochaska, J. X., Bloom, J. S., & 
Thompson, I. B. 2005, ApJ, 634, L25 

Covino, S., Campana, S., Conciatore, M. L., et al. 

2010, A&A, 521, A53+ 

D'Avanzo, P., Perri, M., Fugazza, D., et al. 2010, 
A&A, 522, A20+ 

de Ugarte Postigo, A., Gorosabel, J., Jelinek, M., 
et al. 2006, GRB Coordinates Network, 5290, 1 

Diercks, A. H., Deutsch, E. W., Castander, F. J., 
et al. 1998, ApJ, 503, L105+ 

Douglas, L. S., Bremer, M. N., Lehnert, M. D., 
Stanway, E. R., & Milvang- Jensen, B. 2010, 
MNRAS, 409, 1155 

Eliasdottir, A., Fynbo, J. P. U., Hjorth, J., et al. 
2009, ApJ, 697, 1725 

Evans, P. A. 2011, GRB Coordinates Network, 
12250, 1 

Evans, P. A., & Osborne, J. P. 2011, GRB Coor- 
dinates Network, 12273, 1 

Evans, P. A., Beardmore, A. P., Page, K. L., et al. 
2007, A&A, 469, 379 

— . 2009, MNRAS, 397, 1177 

Ferrero, P., Klose, S., Kann, D. A., et al. 2009, 
A&A, 497, 729 

Finkelstein, S. L., Papovich, C, Salmon, B., et al. 

2011, ApJ, submitted, arXiv:1110.3785 

Fox, A. J., Ledoux, C, Vreeswijk, P. M., Smettc, 
A., & Jaunsen, A. O. 2008, A&A, 491, 189 

Fukugita, M., Shimasaku, K., & Ichikawa, T. 
1995, PASP, 107, 945 

Fynbo, J. P. U., M0ller, P., Thomsen, B., et al. 

2002, A&A, 388, 425 

Fynbo, J. P. U., Jakobsson, P., M0ller, P., et al. 

2003, A&A, 406, L63 

Fynbo, J. P. U., Jakobsson, P., Prochaska, J. X., 
et al. 2009, ApJS, 185, 526 



Fynbo, J. P. U., Laursen, P., Ledoux, C, et al. 
' 2010, MNRAS, 408, 2128 

Fynbo, J. U., Gorosabel, J., Dall, T. H., et al. 
2001, A&A, 373, 796 

Giavalisco, M., Koratkar, A., & Calzetti, D. 1996, 
ApJ, 466, 831 

Gronwall, C, Ciardullo, R., Hickcy, T., et al. 
2007, ApJ, 667, 79 

Grove, L. F., Fynbo, J. P. U., Ledoux, C, et al. 
2009, A&A, 497, 689 

Hamuy, M., Suntzeff, N. B., Heathcote, S. R., 
et al. 1994, PASP, 106, 566 

Hamuy, M., Walker, A. R., Suntzeff, N. B., et al. 
1992, PASP, 104, 533 

Hansen, M., & Oh, S. P. 2006, MNRAS, 367, 979 

Harrington, J. P. 1973, MNRAS, 162, 43 

Hayes, M., Ostlin, G., Mas-Hesse, J. M., et al. 
2005, A&A, 438, 71 

Hjorth, J., & Bloom, J. S. 2012, in Gamma- 
ray Bursts, eds. Kouveliotou, Wijers & 
Woosley, Cambridge University Press, in press, 
arXiv:1104.2274 

Hjorth, J., Sollerman, J., M0ller, P., et al. 2003, 
Nature, 423, 847 

Hjorth, J., Malesani, D., Jakobsson, P., et al. 2012, 
ApJ, submitted, arXiv:1205.3162 

Holland, S. T., Weidinger, M., Fynbo, J. P. U., 
et al. 2003, AJ, 125, 2291 

Jakobsson, P., Hjorth, J., Fynbo, J. P. U., et al. 
2004a, ApJ, 617, L21 

Jakobsson, P., Malesani, D., Hjorth, J., Fynbo, 
J. P. U., & Milvang- Jensen, B. 2011, As- 
tronomische Nachrichten, 332, 276 

Jakobsson, P., Tanvir, N., Jensen, B. L., et al. 
2006a, GRB Coordinates Network, 5298, 1 

Jakobsson, P., Hjorth, J., Ramirez-Ruiz, E., et al. 
2004b, New A, 9, 435 

Jakobsson, P., Hjorth, J., Fynbo, J. P. U., et al. 
2004c, A&A, 427, 785 



24 



Jakobsson, P., Bjornsson, G., Fynbo, J. P. U., 
et al. 2005, MNRAS, 362, 245 

Jakobsson, P., Fynbo, J. P. U., Ledoux, C, et al. 
2006b, A&A, 460, L13 

Jakobsson, P., Hjorth, J., Malesani, D., et al. 2012, 
ApJ, in press, arXiv:1205.3490 

Jaunsen, A. O., Rol, E., Watson, D. J., et al. 2008, 
ApJ, 681, 453 

Kann, D. A., Klose, S., Zhang, B., et al. 2010, 
ApJ, 720, 1513 

Kennicutt, R. C. 1998, ARA&A, 36, 189 

Krimm, H. A., Granot, J., Marshall, F. E., et al. 

2007, ApJ, 665, 554 

Kriihlcr, T., Kiipcii Yoldass, A., Greiner, J., et al. 

2008, ApJ, 685, 376 

Kriihlcr, T., Greiner, J., Schady, P., et al. 2011, 
A&A, 534, A108 

Kriihlcr, T., Malesani, D., Milvang- Jensen, B., 
et al. 2012, ApJ, submitted, arXiv:1205.4036 

Kulkarni, S. R., Djorgovski, S. G., Ramaprakash, 
A. N, et al. 1998, Nature, 393, 35 

Landolt, A. U. 1992, AJ, 104, 340 

Laursen, P., Sommer-Larsen, J., & Andersen, 
A. C. 2009, ApJ, 704, 1640 

MacFadyen, A. L, & Woosley, S. E. 1999, ApJ, 
524, 262 

Malesani, D., et al. 2012, ApJ, in prep. 

Markwardt, C. B. 2009, in Astronomical Society 
of the Pacific Conference Series, Vol. 411, As- 
tronomical Society of the Pacific Conference 
Series, ed. D. A. Bohlender, D. Durand, & 
P. Dowler, 251-254 

McLinden, E. M., Finkelstein, S. L., Rhoads, J. E., 
et al. 2011, ApJ, 730, 136 

Michalowski, M. J., Kamble, A., Hjorth, J., et al. 
2012, ApJ, submitted, arXiv:1205.4239 

Milvang- Jensen, B., Noll, S., Halliday, C, et al. 
2008, A&A, 482, 419 



M0llcr, P., Fynbo, J. P. U., Hjorth, J., et al. 2002, 
A&A, 396, L21 

More, J. J. 1978, in Lecture Notes in Mathematics, 
Vol. 630, Numerical Analysis, Proceedings of 
the Biennial Conference Held atDundec, June 
28-July 1, 1977, ed. G. A. Watson, 105-116 

Neufeld, D. A. 1991, ApJ, 370, L85 

Niino, Y., Totani, T., & Kobayashi, M. A. R. 2009, 
ApJ, 707, 1634 

Nilsson, K. K., Tapken, C, M0ller, P., et al. 2009, 
A&A, 498, 13 

O'Donncll, J. E. 1994, ApJ, 422, 158 

Oke, J. B. 1990, AJ, 99, 1621 

Pandey, S. B., Castro-Tirado, A. J., McBreen, S., 
et al. 2006, A&A, 460, 415 

Patat, F. 2003, A&A, 400, 1183 

Perley, D. A., Cenko, S. B., Bloom, J. S., et al. 
2009, AJ, 138, 1690 

Pettini, M., Stcidel, C. C, Adelberger, K. L., 
Dickinson, M., & Giavalisco, M. 2000, ApJ, 528, 
96 

Piranomonte, S., Ward, P. A., Fiore, F., et al. 

2008, A&A, 492, 775 

Press, W. H, Teukolsky, S. A., Vetterling, W. T., 
& Flannery, B. P. 1992, Numerical recipes in 
FORTRAN. The art of scientific computing 
(Cambridge: Cambridge University Press, 2nd 
ed.) 

Prochaska, J. X., Chen, H.-W., Bloom, J. S., et al. 
2007, ApJS, 168, 231 

Rauch, M., Haehnelt, M., Bunker, A., et al. 2008, 
ApJ, 681, 856 

Rol, E., Jakobsson, P., Tanvir, N., & Levan, A. 
2006, GRB Coordinates Network, 5555, 1 

Salvaterra, R., Delia Valle, M., Campana, S., et al. 

2009, Nature, 461, 1258 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 
1998, ApJ, 500, 525 



25 



Schwarz, H. E., & Melnick, J. 1993, in The ESO 
Users' Manual, 19-24 

Shaplcy, A. E., Steidel, C. C, Pcttini, M., & Adcl- 
berger, K. L. 2003, ApJ, 588, 65 

Stanek, K. Z., Matheson, T., Garnavich, P. M., 
et al. 2003, ApJ, 591, L17 

Stark, D. P., Ellis, R. S., & Ouchi, M. 2011, ApJ, 
728, L2+ 

Steidel, C. C, & Hamilton, D. 1993, AJ, 105, 2017 

Tanvir, N., Rol, E., & Hewett, P. 2006, GRB Co- 
ordinates Network, 5216, 1 

Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, 
Nature, 461, 1254 

Thone, C. C, Perley, D. A., Cooke, J., et al. 2007, 
GRB Coordinates Network, 6741, 1 

Thone, C. C, Kann, D. A., Johannesson, G., et al. 
2010, A&A, 523, A70+ 

Troja, E., Cusumano, G., O'Brien, P. T., et al. 
2007, ApJ, 665, 599 

Tug, H. 1977, The Messenger, 11, 7 

Valls-Gabaud, D. 1993, ApJ, 419, 7 

van Dokkum, P. G. 2001, PASP, 113, 1420 

Verhammc, A., Schaerer, D., & Maselli, A. 2006, 
A&A, 460, 397 

Vreeswijk, P. M., Ellison, S. L., Ledoux, C, et al. 
2004, A&A, 419, 927 

Wainwright, C, Berger, E., & Penprase, B. E. 
2007, ApJ, 657, 367 

Watson, D., Fynbo, J. P. U., Ledoux, C, et al. 
2006, ApJ, 652, 1011 

Wijers, R. A. M. J., Bloom, J. S., Bagla, J. S., & 
Natarajan, P. 1998, MNRAS, 294, L13 

Wooslcy, S. E., & Bloom, J. S. 2006, ARA&A, 44, 
507 

Woosley, S. E., & Heger, A. 2006, ApJ, 637, 914 
Yoon, S., & Langer, N. 2005, A&A, 443, 643 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



26 



A. Observations obtained of systems not in the TOUGH sample 



Three systems which were not in the final TOUGH sample were also observed, see Table [6] The reason 
for these 3 systems not being in the TOUGH sample are as follows: GRB 050603 and GRB 060223A did not 
have an XRT position distributed within 12 hours (although an XRT observation had been made within 12 
hours), and GRB070810A did not have a Sun distan c e grea ter than 55° (its Sun distance was 49°). The 
TOUGH sample criteria are described in iHiorth et al.l (l2012h and are summarised in W2.1\ 

The spectra are shown in Fig. [12] (2D spectra), Fig. [13] (ID spectra) and Fig. [14] (spatial profiles). For 
none of these systems neither the continuum nor the Lya emission line were detected, see Table [7] 

For GRB 050603 the afterglow redshift of z = 2.821 from iBerger fc Becken (|2005l ) is likely wrong: it 
was derived based on a reported bright emission line interpreted as Lya in the afterglow spectrum (0.75 hr 
exposure with Magellan/IMACS), but in our deep host spectrum (2.2 hr exposure with VLT/FORS1) we do 
not detect any emission; we derive a 3cr upper limit on the Lya flux at z = 2.821 of 4.7 x 10~ 18 ergs -1 cm~ 2 
(cf. Table [7]). We do not find an emission line at any other redshift. 
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Fig. 12. — 2D spectra for the 3 systems not in the TOUGH sample. See Fig.UJfor further information. 
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Fig. 13. — ID spectra for the 3 systems not in the TOUGH sample. See Fig. [5] for further information. 
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Table 6 

Observed systems not in the TOUGH sample 



Name 


z 


Ref. 


-Rhost 


Grism+filtcr 


CCD 


fr^total 
exp 


Seeing 


A v 








(mag) 






(hr) 


(arcsec) 


(mag) 


GRB 050603 


N/A a 


(1) 


> 26.6 


600B 


old 


2.2 


< 1.1 


0.092 


GRB 060223A 


4.406 


(2) 


> 26.3 


600R+GG435 


old 


2.1 


0.7 


0.389 


GRB 070810A 


2.17 


(3) 


> 26.7 


600B 


new 


1.5 


0.7 


0.072 



a The redshift z = 2.821 for GRB 050603 from iBerger fc Beckerl |2005l ) is likely wrong: it was 
derived based on a reported very bright emission line interpreted as Lya in the afterglow spectrum, 
but in our deep host spectrum we do not detect any Lya emission; we derive a 3<r upper limit on 
the Lya flux at z = 2.821 of 4.7 x 10~ 18 ergs" 1 cm" 2 (cf. Table[7|). 



Note. — See Table[T]for further information. 

References. — m lBerger fc Becked (|2005h ; (21 ICharv et~ai1 |2007l ); f3) lThone etHI (j2007T ). 



Table 7 

LYa MEASUREMENTS from the spectra 



Name 




Lya aperture 




F(Lya) 


L(Lya) 


F A (cont.) 


EW(Lya) 


v(Lya) 




c(v) 


c(s) 




w(s) 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


GRB 050603 


281 


0.02 


906 


1.20 


< 9.3 


< 0.62 


< 9.0 






GRB 060223A 


297 


0.00 


862 


1.20 


< 14.5 


< 2.80 


< 12.6 






GRB 070810A 


264 


-0.05 


912 


1.25 


< 8.7 


< 0.31 


< 9.5 







Note. — See Table [3] for further information. 
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Fig. 14. — Spatial profiles of the 3 systems not in the TOUGH sample. See Fig. [6] for further information. 
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